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A combined-isotopologue direct-potential-fit (DPF) analysis of optical and pho- 
toassociation spectroscopy data for the a 3 X+ and 1 3 S+ states of Li2 has 
yielded accurate analytic potential energy functions for both states. The recom- 
mended M3LRg'3(3) potential for the a 3 S+ of 7,7 Li2 has a well depth of D e = 
333.758(7) cm -1 and equilibrium distance of r e = 4.17005(3) A, and the associated 
scattering lengths are agh = —14.759(9) A for 7,7 Li2 and —1906(50) A for 6 ' 6 Li2. For 
the 1 3 S+ state, in spite of a gap of ~ 5200 cm^ 1 (from w(l 3 S+) = 8 — 61) for 
which there are no data, the DPF procedure has no difficulty determining an accu- 
rate overall potential. The 1 3 E+ state of the ' 7 Li2 isotopologue has a well depth of 
D e = 7092.417(33) cm" 1 and equilibrium distance of r e = 3.06524(9) A. The long- 
range tail of the recommended M3LRgg(9) potential energy function for the l 3 Sj" 
state is defined by the lowest eigenvalue of a 3 x 3 long-range interstate coupling 
matrix to take into account the 3-state mixing near its asymptote. 
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I. INTRODUCTION 



Modern theoretical studies of ultra-cold atomic gases demand a very accurate knowledge of the 
potential energy curves (PECs) of the systems of interest. Since Li 2 is the second smallest uncharged 
stable homonuclear molecule, its chemical and physical properties are particularly interesting. In 
recent years, considerable effort has been focussed on the lowest singlet states of Li 2 .-~— However, 
the properties of the triplet states of Li 2 are much less well known. 

The first observation of discrete spectra involving the lowest triplet state of Li 2 was reported 
in 1985 by Xie and Field,— who used perturbation-facilitated optical-optical double resonance 
(PFOODR) techniques to excite 2 3 IT g — a 3 ^u emission. They observed transitions involving 
f(a 3 S+) = — 6, but because of the limited resolution available at the time, their results have 
been superceded by later work. The first high-resolution triplet-system measurements were reported 
three years later by Martin et al.,— who performed a Fourier transform study of the 1 3 £+ — > a 3 E+ 
system of 7 ' 7 Li 2 involving v' = 1 — 7 of the upper state and v " = — 7 of the ground triplet state, 
with average measurement uncertainties of only ±0.01 cm± Analogous results for the same system 
of 6 ' 6 Li 2 , spanning the same ranges of vibrational levels, and with the same accuracy, were reported 
a year later by Linton et al. 13 A decade later Linton et alr^ reported a high-resolution version of 
the PFOODR experiment of Xie and Field^ which yielded accurate (uncertainties ranging from 
0.005 to 0.001 cm -1 ) observations of transitions into u(a 3 E+)=0 — 9 from a handful of rotational 
levels of the v' = 1 and 2 levels of the 2 3 Il g state of 7 ' 7 Li 2 . In addition , a two-photon photoassoci- 
ation spectroscopy (PAS) experiment by Abraham et al— had yielded a direct measurement of the 
0.416(±0.001) cm' 1 binding energy of the v=10, N=0 level of the a 3 £+ state of 7 ' 7 Li 2 . 

In addition to the seven vibrational levels v = 1—7 of the 1 3 £+ state observed in the emission 
experiments, the binding energies of levels v = 62 — 90 of 7 ' 7 Li 2 and v = 56 — 84 of 6 ' 6 Li 2 were 
measured in a PAS study by Abraham et al— However, to date there has been no reported attempt 
to bridge the 5100 cm -1 gap between the two sets of results in order to provide a global description 
of this state. This problem is illustrated by Fig. 1, which shows the regions of the a 3 S+ and 1 3 £+ 
potentials associated with the currently available experimental data. The task of bridging the chasm 
for the 1 3 £+ state is complicated by the fact that in the region very near the upper-state asymptote 
spanned by the PAS data, a transition from case (b) to case (c) coupling leads to a mixing of the 
1 3 E+ state with two other states that also have l g symmetry in this long-range (r > 20 A) region. 

The only PECs which have been reported for the lowest triplet states (a 3 S+ and 1 3 S+) of Li 2 , 
are point-wise, semiclassical RKR curves generated from Dunham or near-dissociation expansions 
for the vibrational energies and inertial rotation B v constants.—"— For the 1 3 E+ upper state, those 
PECs were based only on data for vibrational levels v ' = — 7 ,— ^ since that work preceded the the 
photo-association spectroscopy (PAS) studies of this system.—^ Thus, the best available potential 
for the 1 3 E+ state provides no realistic predictions for those subsequently observed high vibrational 
levels, took no account of interactions with other states of Li 2 near the dissociation asymptote 
of l 3 Ejj~, and did not incorporate the theoretically known inverse-power long-range behaviour. 
Finally, in all studies of these states to date, 7 ' 7 Li 2 and 6,6 Li 2 were treated independently, and as a 
consequence, the effect of Born-Oppenheimer breakdown (BOB) in this system remains unknown. 
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The present work presents a fully quantum mechanical direct-potential-fit (DPF) data analysis 
which accounts for all of the optical and PAS data described above in terms of global analytic 
potential energy functions for the ot 3 S+ and 1 3 £+ states of Li 2 , while independent term values are 
used to represent the 20 levels of the 2 3 IT 5 state giving rise to the observed high-resolution emission 
into v(a 3 E+) = — 9. The three longest-range inverse-power contributions to the interaction 
energy are incorporated into the a 3 £+ and 1 3 Ejj" state potential energy functions, and the function 
for the 1 3 E+ state explicitly accounts for the three-state mixing referred to above. In addition, 
the incorporation of adiabatic Born-Oppenheimer breakdown (BOB) correction functions in the 
Hamiltonians for the a 3 £+ and 1 3 S+ states allows the data for the two isotopologues to be treated 
simultaneously. 

Three aspects of this system made its analysis unusually challenging. Firstly, the long-range tail 
of the 1 3 £+ state PEC is not the familiar sum of simple inverse-power terms, since the 1 3 S+ state 
couples strongly to two other states near the dissociation asymptote (see § II.C.l). Secondly, the 
fact that the leading long-range term of the PEC tail for the 1 3 E+ state is C3/V 3 means that the 
Morse /long-range (MLR) potential function leads to unphysical long-range behavior if not addressed 
appropriately (see § II. C. 3). 10 Finally, as illustrated by Fig. 1, the data for the 1 3 E+ state has gaps 
from v = 7 to v = 62 for 7 ' 7 Li 2 and from v — 7 to v — 56 for 6,6 Li 2 , which span more than 72% of 
the well depth. Such a large gap in experimental information has never (to our knowledge) been 
treated successfully by a potential-fit analysis in a purely empirical manner. 



II. MODELS AND METHODOLOGY 

A. DPF Data Analyses and the Radial Hamiltonian 

In a DPF spectroscopic data analysis, the upper and lower level of every observed energy tran- 
sition is assumed to be an eigenvalue of an effective radial Schrodinger equation characterized by 
a parameterized potential energy function and (when appropriate) parameterized radial strength 
functions characterizing appropriate BOB terms. Given some plausible initial trial parameter val- 
ues for characterizing the relevant potential, solution of the associated Schrodinger equation yields 
an eigenvalue E Vi j and eigenfunction i/j v ,j(r) for each observed level. The difference between the 
energies of appropriate upper and lower levels then yields an estimate of each observed transition 
energy, while use of the Hellmann-Feynman theorem: 



dE Vy 
d Pj 



dH 

dpi 



(1) 



yields the partial derivatives required for performing a least-squares fit of the simulated transitions 
to the experimental data. 

Since the observed transition energies are not linear functions of the parameters defining the 
effective radial Hamiltonian, a DPF analysis requires the use of an iterative non-linear least-squares 
fitting procedure. The quality of a given fit is characterized by the value of the dimensionless 
root-mean-square deviation of the N experimental data y bs(i) from the predicted values y ca Xc(i) 
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generated from the relevant Hamiltonian(s): 
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in which u(i) is the uncertainty in the reported value of experimental datum i. In the present work, 
these fits were performed using the publicly available program DPotFit, 17 while the requisite 
initial estimates of the potential function parameters were obtained using the program betaFIT 18 
with preliminary RKR potentials generated using conventional Dunham expansions. 

As with most diatomic DPF analyses reported to date, the present work is based on the effective 
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radial Schrodinger equation presented by Watson,—^ and uses the conventions described in Refs. 
and |22j. In particular, the rovibrational levels of isotopologue a of diatomic molecule AB in a given 
electronic state are the eigenvalues of the radial Schrodinger equation: 
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in which (r) is the effective adiabatic internuclear potential for a selected reference isotopologue 
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is the difference between the effective adiabatic 
r) is the non-adiabatic centrifugal-potential 
correction function for isotopologue-a, and the reduced mass \i a is defined by the atomic masses 



potentials for isotopologue-a and isotopologue-1, g( a > 



M^ and M B . Following standard conventions,— — the BOB terms Al^ (r) and g^ a \r) are both 
written as a sum of contributions from component atoms A and B: 
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Mjy B — M^j B are the differences between the atomic masses of atoms A or 



B in isotopologue-a and in isotopologue-1. In the present case A = 
collapse to 



B = Li, and these expressions 
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aA yr) and i?^(r)), two mass factors must be retained in order to allow us to describe all possible 



Although only a single radial strength function of each type must be considered in the present case 
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molecular isotopologues. 
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B. The 'Basic' Morse /Long- Range (MLR) Potential Energy Function 

The next step is to introduce an optimal analytic function for representing the effective adiabatic 
internuclear potential for the reference isotopologue, V^\r) = V(r) . The present work is based 



on use of the version of the Morse/long-range (MLR) potential energy function of Refs.llOl and 
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VmLrM = D e 



(8) 



in which T) e is the well depth, r e the equilibrium internuclear distance, and the radial variable in 
the exponent is 
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The parameterized exponent coefficient function /3(r) which governs the details of the shape of the 
potential is defined so that 

lim 0{r) = /3oo = In ( -^-] , (10) 
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and as a result, the long-range behavior of the potential energy function is defined by the function 
«lr(?") : 

V ULR (r) ~ D e - u LR (r) + O (n LR (r) 2 /4D e ) , . . . , (11) 

while the denominator factor u\, R {r e ) is simply the value of that long-range tail function evaluated 
at the equilibrium bond length. 

The theory of long-range intermolecular forces shows us that in general, mlr(^) may be written 
in the form 

uUr) = £ D m (r) ^ • (12) 

in which the powers and coefficients C mi of the terms contributing to this sum are determined by 
the symmetry of the given electronic state and the nature of the atoms to which it dissociates,— _ — 
and the 'damping functions' D m (r) were introduced to take account of the weakening of the inter- 
action energies associated with these simple inverse-power terms due to overlap of the electronic 
wavefunctions on the interacting atoms.~ While most previous applications of the MLR potential 



function form omitted the D m (r) damping function factors, it was shown in Ref.l24i that in addition 
to providing a more realistic physical description of the long-range potential tail, their introduction 
improves the extrapolation behaviour of the repulsive short-range potential wall, and when they are 
included, fewer parameters are required to achieve a given quality of fit to experimental data. In 
either case, the structure of Eq. ffTT]) means that at large distances where wlr(^) 2> «LR(r) 2 /(4 D e ), 
the long-range behaviour of Vmlr is defined by mlr(^). 

Following Ref.|24j, the present work uses the modified Douketis-type^ damping function form 



b ds (s)(pr) - ds '- '-- 2 \ m+1 



D™( s \ r ) =11- e « , (13) 



with s = — 1. Here, b ds (s) and c ds (s) are treated as system-independent parameters with 6 ds (— 1) = 
3.30 and c ds (— 1) = 0.423 . For interacting atoms A and B, p = pab = 2/?aPb/(pa + Pb), in which 
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Pa = (Ip I 'I^) 2 ^ 3 is denned in terms the ionization potential of the atom A and that of an H 
atom (Ip and 1^ respectively). Inclusion of these damping functions means that at very short 
distances Vmlr^) oc 1/r 2 .— Comparisons with ab initio results for a sampling of chemical and 
Van der Waals interactions showed that this type of damping function yielded quite realistic MLR 
short-range extrapolation behaviour,™ so this D m (r) form is adopted here. 

In order to ensure that the exponent coefficient function in Eq. (jSJ) satisfies Eq. (flQ|) . it is cus- 
tomary to write it as the constrained polynomial: 
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0(r) = 0}(r) = yf{r)^ + [l - yf (r)] J> (r)]* . (14) 

8=0 

This function is expressed in terms of two radial variables which are similar to yp q {r), but are defined 
with respect to a different expansion center (rw), and involve two different powers, p and q (the 
reasons for this structure are discussed in Ref.llOl): 



<M = -y— -^7 and yf (r) ee . (15) 

The limiting long-range behaviour of the exponential term in Eq. (JH]) gives rise to additional inverse- 
power contributions to the long-range potential of Eq. (11 ip . with the leading term being proportional 
to \/r mi+p . This means that the power p must be greater than (mi ast — mi) if the long-range 
behaviour of Eq. (I12p is to be maintained.— There is no analogous formal constraint on the value 
of q ; however, experience suggests that its value should lie in the range 2 < q < p In 
early work with this potential function form, the radial variables in Eq. ( I14p were both defined as 
yp^ir) of Eq. fl9]) (i.e., r ref = r e and q = p)M>^— However, it has since been shown that setting 
r ref > r e and q < p can significantly reduce the number of parameters required to describe a 
given data set accurately, and yields more stable expansions.— >2a22 

A second consideration associated with the use of the damping functions of Eq. ([13]) is their effect 
on the shape of the short-range repulsive potential wall of an MLR potential. As was pointed out 



in Ref.l24l the fact that the radial variables ylfJr) — > — 1 as r — > means that at very small 



distances Vmlr(^) oc {mlr(^)} 2 - If damping functions are neglected (i.e., if Eq. (fT2|) did not include 
the D m (r) functions), then the limiting short-range behaviour of the potential energy function would 
be Vmlr(^) oc l/r 2miast . However, for a typical two- or three-term ^lr(^) expansion, mi ast = 8 
or 10, and the resulting 1/r 16 or 1/r 20 short-range repulsive wall behaviour would be unphysically 
excessively steep. In the data-sensitive region of the potential well, this excessive growth rate would 
be compensated for by the behaviour of the empirically determined exponent coefficient function 
f3(r). However, the unphysical high-order r -16 or r -20 singular behaviour would re-assert itself in 
the shorter-range extrapolation region. 

In this paper, the label for particular MLR potential function models is written as MxLRp r ° f (N), 
in which x is the number of inverse-power terms incorporated into Ulr(^), while p, q, r re f and 
are defined above. The 'basic' MLR model described above is used herein to describe the potential 
energy function for the o 3 S+ state of Li 2 . However, some enhancements were required for treating 
the 1 3 S+ state. 



C. Modified MLR Potential for the 1 3 S+ State of Li 2 

1. Incorporating interstate coupling into the MLR model 

At a preliminary stage of the present work, the 1 3 E+ state of Li2 was represented by the 'basic' 
MLR function of Eqs. fl8l- ffl~5l) . in which u\^{r) consisted of three terms, with rrii = {3,6,8}. 
Because of the added complexity due to interstate coupling near the dissociation asymptote, all 
damping functions for this state were fixed at D m (r) = 1. This model was able to provide an 
excellent fit both to the fluorescence data for v' = — 7, and to the 7 ' 7 Li2 PAS data for v = 62 — 70 
and the 6 ' 6 Li2 PAS data for v = 56 — 65 whose upper limits which correspond to binding energies of 
about 24 cm -1 . However, when PAS data for higher vibrational levels were included in the analysis, 
the quality of fit got progressively worse, and the discrepancies could not easily be removed simply 
by increasing the order of the the polynomial /3(r). The reason for this increasing inability of the 
basic MLR model to account for levels lying very near dissociation is that the 1 3 £+ state of Li 2 
couples to two other states near its dissociation asymptote. 

This same type of problem was encountered in a recent study of the A( 1 E+) — X( 1 S+) system of 
Li 2 . In that case the 0+(A 1 E+) state which goes to the Li( 2 Py 2 ) + Li( 2 Si/ 2 ) asymptote couples to 
the 0+(fe 3 IT) state which goes to the higher Li( 2 Py 2 ) + Li( 2 S'i/ 2 ) limit,— ^ and the energies of levels 
lying near dissociation could not be explained properly without taking account of the inter-state 
mixing. Fortunately, Aubert-Frecon and co-workers had derived an analytic description of those 
coupled states based on the eigenvalues of a 2 x 2 interaction matrix,-^ and it was shown that their 
analytic expression for the lower eigenvalue could readily be used to define «lr(^) for this state in 
an MLR potential model.— 

Treatment of levels lying near the dissociation limit of the 1 3 £+ state of Li 2 involves a similar 
problem; while it dissociates to the Li( 2 Pi/ 2 ) + Li( 2 5'i/ 2 ) limit, it couples to l 9 ( 1 n) and l g ( 3 II) states 
which correlate with the higher Li( 2 Py 2 ) + Li ( 2 Si /2 ) limit.- 3 ^ Since the Li( 2 P) state spin-orbit 
splitting is quite small (ca. 0.335 cm -1 ), the interstate coupling only becomes important for levels 
lying relatively close to the dissociation limit. Fortunately, Aubert-Frecon and co-workers have 
studied this case too.— In particular, they presented expressions for the six independent elements 
of the symmetric 3x3 matrix that defines the long-range interaction energies for these three states. 
Their matrix elements took into account the first-order resonance-dipole (1/r 3 ) term, the leading 
dispersion energy terms, and the exchange energy. If we neglect the exchange term, keep only the 
first two (m = 6 and 8) dispersion energy terms, set the zero of energy at the 1 3 E+-state asymptote, 
make use of the symmetry relation for m = 3 , 

C 3 S ° = 2Cl U9 = -2C3 119 = Cf , (16) 



and that for m — 6 , 

C,^ = CT* = C 6 n , (17) 
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3 S + 

and define Cg 8 = C 6 8 9 , their 3x3 long-range interaction matrix Mlr becomes 

C 3 I Cf+2C^ . Cf + C s 3 +C S 3 \ y/2 I Cf , C^-C" . 2Cf — C a 3 — C s 3 \ _\_ ( C3. 

T 3 " 1 7 s 1 ?S I 3 I r 5 " 1 ^ 1 I 
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\~f^ 1 F» y 2V3 V 7^ 1 r 8 ^ y ^^r -r ^ 2rS - J^ AE J 

(18) 

in which A.E is the accurately known (positive) spin-orbit splitting energy of Li( 2 P).— Note that 
in contrast with Ref.l37l. the present formulation treats attractive C m coefficients as positive, rather 
than negative quantities. Following the correlation scheme given by Movre and Pichler,— the lowest 
eigenvalue of the matrix (fT8l) defines the long-range tail of the 1 3 E+-state interaction potential; see 
Fig. 2. 



Analytic expressions for the three eigenvalues of Eq. ( fT8l) were reported in Ref.|37J. Those ex- 
pressions for the zeros of the characteristic polynomial for Mlr were obtained using the method 
of Scipione del Ferro and Niccolo Fontana Tartaglia (first published 39 by Gerolamo Cardano in 
1545), and by applying a trigonometric substitution to avoid expressions involving square- roots of 
negative quantities. However, Kopp has demonstrated that while it is useful for obtaining symbolic 
expressions, this formula can yield substantial errors when used for actual computations, primarily 
because of the numerical errors that accumulate when computing the arctan function within the 
formula.— Moreover, the symbolic expressions for the derivatives of the lowest eigenvalue with 
respect to the Cn^ U coefficients required by the least-squares fitting procedure are inconveniently 
complex. Because of these problems, in the present work the eigenvalues of this interaction matrix 
were calculated numerically (using the Jacobi eigenvalue algorithm^^i), and their derivatives with 
respect to the C n coefficients were computed using the discrete version of the Hellmann-Feynman 
theorem: 



dp 



dM 
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dp 

in which Aj is the appropriate eigenvalue of the matrix Mlr, and <j)\ i is the corresponding eigen 



vector. 



2. Simplifying the treatment of interstate coupling for Xi2(l S 



The treatment of the long-range behaviour of the A 1 ^ state of LI2 in Ref.llOl was precisely 
analogous to that for the 1 3 S+ state discussed here, except that while that case involved a 2x2 
matrix whose eigenvalues were determined analytically, the present case involves the 3x3 matrix 
of Eq. (fT8j) whose eigenvalues are calculated numerically. In the treatment of the A x £+ state, 
it was shown that the C^/r 6 and C^/r 8 terms had virtually no effect on the lower (S-state) 
eigenvalue of the 2x2 long-range interstate coupling matrix. This led us to consider making the 
same simplification here. 



Following the approach of Ref . (jlOl ) , we compared the values of the lowest eigenvalue of Eq. (pi 
obtained when all C^ U and coefficients were defined by the theoretical values of Tang et 

al.,— with those obtained on setting = Cg 1 * 9 = C 8 Ub = . Over the range r = 2 to 500 A, the 
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difference between these two estimates of the lowest eigenvalue were always less than 3 x 10 6 cm 1 . 
Thus, it seems clear that in the present treatment of the 1 3 E+ state of Li 2 , no significant errors will 
be introduced if contributions involving Cg 1 , C 8 ns and C 8 Ug are omitted from Eq. (lisp . At the same 
time, it is important to note that these and Cg 1 coefficients cannot be neglected when using the 
two higher eigenvalues of Eq. ( I18p to define the long-range tails of the 3 fl 9 and 1 H g states which 
couple with the 1 3 X+ state of interest here. This point is illustrated by Fig. 3, which compares plots 
of the three eigenvalues of Eq. ( 1181) obtained using all of the Cn^ U coefficients of Tang et ai— (solid 
red curves) with those obtained from this same matrix when = = Cg 1 * 9 = C 8 Ub (dashed blue 
curves). It is clear that at the smaller distances where the Cq and Cg terms become important, one 
cannot use the above approximation when calculating the eigenvalues of ( ITS]) associated with the 
two Ug states. 



3. Implications of the quadratic term in the MLR potential function form 



It was shown in Ref.|l0| that contributions from the quadratic term in Eq. ( II ip can give rise 



to spurious perturbations in long-range behavior of the MLR potential function form. In the 
present case, the leading terms in the long-range potential for the 1 3 S+ state of Li 2 correspond to 
m, = {3,6,8}. If we temporarily ignore the effects of damping and interstate coupling in order to 
write mlr(^) as a simple inverse-power sum, the presence of the quadratic term in Eqs. flS} and (TIT]) 
mean that the effective long-range behavior of the MLR potential would be 

iy\J ^"U iyC) ty\J lyU 

Thus, if the overall effective long-range behavior is to be defined by an inverse-power sum governed 
by the specified C 3 , C 6 and C 8 coefficients (and not include the last two terms in Eq. (120]) !). the 
definition of mlr(^) must compensate for the quadratic terms by being defined as 
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in which Cf =C e + (C 3 ) 2 /(4D e ) and Cf = C 3 Cf/(2®,). 

Since the long-range tail of our 1 3 S+-state potential also includes interstate coupling, these 
expressions for and (which were derived analytically for potentials with simpler long- 
range tails) need to be tested in order to determine how well they cancel the effect of the spurious 
last-two terms in Eq. ( 120 p . Results of a numerical test of this question are presented in Fig. 4, which 
displays plots of the quantity C| (r) = r 3 [D — mlr(^)] vs. r~ 3 for three different definitions of 
wlrO")- This type of plot illustrates the nature of the long-range interaction on a reduced scale. 
If «lr(^) was defined as the simple inverse-power sum of Eq. (121) . as r~ 3 — > the resulting plot 
would approach an intercept of C 3 with a limiting slope of C| ff . In the present case, however, the 
three-state coupling near the potential asymptote causes all of the plots in Fig. 4 to drop off sharply 
for r- 3 < 10~ 4 A~ 3 . 

The desired long-range behavior is achieved when wlrM is defined simply as the lowest eigen- 
values of the matrix of Eq. (ITS]) , mlr(^) = — A m j n (r). For the case in which = C 8 n9 = C 8 n9 = 
and the other Cn^ U coefficients are fixed at the values of Tang et al., 42 this desired behavior is 
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defined by the solid black curve in Fig. 4. The dash-dot-dot red curve in Fig. 4 then shows how the 
long-range behavior of the associated MLR potential, which includes the quadratic term of Eq. (11 ip . 
deviates from this desired long-range behavior. Next, the dotted blue curve shows the effect on the 
long-range MLR behavior of replacing by the quantity defined above. It is immediately 
clear that this removes most of the discrepancy with the 'ideal' long-range behavior (solid black 
curve). Finally, the dashed green curve shows the effect on the long-range MLR potential tail of 
also including the Cg^/r 9 term in the definition of mlr(^) in order to cancel out the spurious r~ 9 
term in Eq. (1211) . It is clear that use of the resulting definition 

u LR (r) = - X min (C 3 ,Cf,C 8] r) + Cf/r 9 (22) 

brings the long-range tail of the overall MLR potential function into essentially exact agreement 
with the desired form. 



4- Inclusion of retardation in the model potential for Li2(l 3 Sjj~) 

It has long been known that at the very large distances where the C^/r 3 term comes to dominate 
the interaction energy in this type of system, "retardation" effects due to the finite speed of light 
should not be neglected.— 1 ^ It was shown by Meath^ that the effect of retardation on an s/p 
resonance-dipole interaction can be accounted for by multiplying Cf by the function f^ t (r) and Cg 1 
by the function f^ t (r), where 
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(23) 
(24) 



in which \ S p — ^sp/^ 71 an d Xsp is the wavelength of light associated with the atomic 2 S — 2 P 
transition. 

It is a straightforward matter to incorporate this retardation behavior into the MLR potential 
function form. In particular, on setting C^ 1 = = 0, making use of the symmetry relationships 
among the Cf, C^ 3 and C^ 3 coefficients, and replacing C$ by C§ )&A \ the long-range interstate 
coupling matrix for the three l g states dissociating to yield Li( 2 Sy 2 ) + Li( 2 P) becomes 
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(25) 



Unless stated otherwise, the definition of the long-range tail of the MLR potential for the 1 3 £+ 
used throughout the rest of this study is therefore given by 



u 



ret 
LR 



r) 



ad J /V.9 



(26) 



in which X r ^\ n (Cf, Cg dj , Cf; r) is the lowest eigenvalue of the interaction energy matrix of Eq. (125]) . 
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D. Born-Oppenheimer Breakdown Functions 

The radial strength functions in Eqs. (E]) and (J7|) may be written as polynomials constrained to 
have specified asymptotic values using the format of Eq. (JTj 



S%(r) = yZ(r) t£ + [1 - yZ(r)} £ vZ( r Y ( 27 ) 

i=0 

B£(r) = yZ(r) C + [1 - (r)] £ !^(r) 4 (28) 

i=0 

in which and are the values of these functions in the limit r — > oo , Mg 1 and tg 1 define their 
values at r = r e , and the radial variables are versions of Eq. (JUJ) associated with chosen values of 
the integers p a d, p n a,, <?ad and g na . 22 The discussion of Ref.H shows that = 0.0 for any molecule 
which dissociates to yield an uncharged atom-A, so t^(a 3 £+) = t^(l 3 £+) = 0.0. In addition, 
we adopt the Watson convention of setting the parameter t^ 1 = 0.0 for both the a 3 S+ and 1 3 £+ 
states, since its value cannot be determined from transition-frequency data alone.— £022. 

Following standard conventions,— the absolute zero of energy is defined as the energy of ground- 
state atoms separated at r — > oo , so by definition w^(a 3 £+) = 0.0. Since the 1 3 S+ state of Li 2 



dissociates to one ground-state ( 2 5i/2) and one excited-state ( 2 -Pi/2) atom, the value of -u^(l 3 £+) 
is then defined by the difference between the atomic 2 -Pi/ 2 2 5*1/2 excitation energies for 6 Li and 
7 Li, which is^ 3 ^ 

SE 6 7 H( 2 P 1/2 ) = AE 6L \ 2 P 1/2 ^ 2 S 1/2 ) - AE 7u ( 2 P 1/2 <- 2 S 1/2 ) = - 0.351338 [cm" 1 ] (29) 

This is the difference between the energy asymptotes of the 1 3 E+ states of 6 ' 6 Li2 and 7,7 Li 2 , and it 
defines the asymptotic value of the adiabatic radial strength function.— Since we select 7 ' 7 Li 2 as 
the reference isotopologue, this yields 

ttS(l s S+) = 5l££CP 1/2 ) 1 2 (l - = 1.05574 [cm" 1 ] (30) 

We now address the choice of powers p a d, q&d, Pna and q na for defining the radial expansion 
variables in Eqs. (12 7p and ( 128]) . As was pointed out in Ref.|2l|, if the effective adiabatic potential for 
the 'minor' isotopologues is to have the same limiting long-range behavior as that for the reference 
isotopologue, p a d must be greater than or equal to the power of the longest-range term in the 
intermolecular potential for that state. Thus, we set p a a(a 3 S+) = 6 and p a d(l 3 S+) = 3. Note that 
BOB radial strength functions are relatively weak and slowly varying, and few terms are required 
to define them. As a result, there is no need here to introduce an r re f ^ r e extension into the 
definition of the expansion variables in Eqs. f l27l) and (j28~l) . and for the sake of simplicity we set 
g ad (a 3 S+) =p ad (a 3 S+) =6 and g ad (l 3 £+) = Pad (l 3 S+) = 3 . 

We are not aware of any theoretical predictions regarding the limiting long-range behavior of the 
centrifugal non-adiabatic radial strength function i?^ d (r), so we have no basis for assigning particular 
values to p na . Moreover, as in the discussion of_§ III.B.(iii), there are no physical constraints on the 



values of g na . At the same time, Fig. 3 of Ref.l22l shows that use of too small values for these powers 



can give rise to physically implausible extrema in the resulting functions on the interval between 
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the data region and the asymptote, while use of too high values will lead to a requirement for an 
excessive number of expansion coefficients. For simplicity, we therefore chose to set p nSL = g na = 3 
in fits to models which included non-zero -R„a( r ) functions for either electronic state. 

Since y^(r e ) = and y V p C f d {r — > oo) = 1, the algebraic form of Eq. ff27|) means that the difference 
between the well depths of different Li2 isotopologues in a given electronic state is given by the 
expression^ 



r ;a I ;b 



and that the analogous shift in the equilibrium distance r e is 

K '' ~ [m$ m$) I ' (32) 

in which k is the harmonic force constant at the potential minimum in units cm -1 A , and 

,,, | ^gad \ (Mqq - u )p ad + Ui g ad , s 

/ r=r e 

Similarly, the electronic isotope shift will be 

5 { AT^ } = AT^ - AT^ (34) 



AMff 



Note that in the present context, D i and are the values of the well depth and equilibrium 
distance of the MLR potential for the reference isotopologue species, and are determined by the 
DPF analysis. 

Finally, as was pointed out by McAlexander et al.,— for the 1 state of LI2, the dominant 
BOB contribution to the rotationless potential at large r has the form 

AV^ * 2B^(r) = 2(^), (35) 

and since \x a for isotopic Lip is relatively small, this behavior must be considered. Following the 



approach of Refs. 
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45 



and 



4(1 we have chosen to treat this term as a separate additive contribution 



to the effective interaction potential for each isotopologue, which therefore takes on the form: 

V^Jr) = V MLR (r) + Atf(r) + 2B«(r) . (36) 

As was pointed out by Vogt et al.,— this AV a d{r) term is readily incorporated into the Hamiltonian 
by simply replacing the factor [J( J+l)] in Eq. fl3]) by [J(J+ 1) +2], and their approach was adopted 
here. However, this means that the overall 1 3 E+-state well depth and equilibrium distance are 
actually D* ot (c) = - 2B^(r e ) and r* ot (c) = r ( e a) + AB^(r e )/(kr e ) , where D ( e a) and r { e a) 
are defined by Eqs. ( )3Tj) and ( )32|) . and 3)i and are the (fitted) reference- isotopologue MLR 
parameters for that state. 
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III. POTENTIALS FOR THE a 3 S+ AND 1 3 S+ STATES OF Li 2 

A. Data Set and Methodology 

An overview of the experimental data used in this work is presented in Table I. Most of the 
data (2555 out of 2792 observations), came from fluorescence experiments performed by Martin et 
alM and Linton et a/.— That data set was enlarged by inclusion of 137 2 3 n s — a 3 £+ transitions 
of 7,7 Li 2 taken from the study of Ref. 1J, which extended the a 3 S+ vibrational range to v = 9. 
Finally, information about levels lying very near dissociation is provided by the one available PAS 
datum for the v = 10 level of the a 3 E+ state,— and 99 PAS data for the higher levels of the 1 3 £+ 
state.— ^ Throughout this study, the upper levels of all transitions originating in the the 2 3 I1 9 
state were represented by independent fitted term values. A listing of the experimental data used in 
the present analysis in included in the Supplementary Data supplied to the journal's www archive 

All of the DPF data-analysis fits described herein were performed using the program DPotFit, 
which is freely available (with a manual) for download from the www.— The initial trial values of the 
parameters /3j required for those fits were generated by applying the program betaFIT (also available 
from the www) 18 to sets of turning points obtained from preliminary versions of the analysis. 

B. Model for Li 2 (a 3 £+) 

The a 3 S+ state of Li 2 dissociates to yield two S'-state atoms, and ignoring hyperfine effects, 
there is no noteworthy interstate coupling. The theory of intermolecular forces therefore tells us 
that the leading contributions to the long-range intermolecular potential should consist of terms 
associated with (inverse) powers = {6,8,10}. The present analysis therefore represented the 
potential energy for this species by an MLR potential incorporating the long-range tail function 

,&*>M = D 6 (r) £ + D sW £ + D m{ r) , (37) 

in which D m {r) are the modified Douketis-type damping functions of Eq. (fl3|) . p = 0.54,~ and the 
dispersion energy coefficients for this state were fixed at the values reported by Tang et ala- 
lia the initial work to determine an optimum model for this state, the PAS data for the 1 3 S+ 
state were ignored and all of the observed levels of both the 1 3 S+ and 2 3 IT 9 states were represented 
by independent term values, so only a single potential energy function was involved in the analysis. 
Fits were then performed to a wide variety of models corresponding to different choices for the order 
N of the polynomial in Eq. (fl4|) . and for the power q and the reference distance r re f of Eq. (fIBl) . As 
was pointed out in § II.B, the power p must be larger than the difference between the largest and 
smallest powers of the terms contributing to Eq. ( 137ft . so it was fixed as p = 5. 

Figure 5 summarizes results obtained for six families of potential function models of this type. 
For any given values of q and r re f, the fits will always converge when the polynomial order N 
becomes sufficiently large. This point is illustrated by the convergence of the four families of solid 
triangular points shown in the lower panel of Fig. 5, which correspond to models with q = 3 and 
N increasing from two to five. The results for N = 3 with q = 2 and 4 (open square and round 
points, respectively) show that the optimum value of r Te{ will depend on other features of the model, 
but convergence to the same limiting quality-of-fit dd and essentially the same limiting values of 
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physically interesting parameters such as D e will be achieved in any case. The upper panel of Fig. 5 
shows the fitted values of D e for the various models considered in the lower panel. Those results 
show that the fitted value of D e may vary considerably from one model to another, but at the r re f 
value where dd approaches its minimum value, all models converge to essentially the same value of 

All of the fits that are summarized by Fig. 5 used all of the 2693 'optical' data for the two 
isotopologues 7,7 Li 2 and 6 ' 7 Li 2 , plus the one PAS datum which had been reported for the a 3 £+ 
state, while the 395 term values of the observed levels of the 1 3 £+ state and 20 term values for the 
2 3 Il g state were treated as free parameters. As might be expected when dealing with light atoms 
such as Li, BOB effects are not negligible in this system. In particular, it was found that allowing 
for one non-zero term u^ 1 in the expression for the 'adiabatic' correction radial strength function of 
Eq. (T2T|) reduced the value of dd by 2.4%; however, freeing a second coefficient (uY) only reduced 
dd by an additional 0.08%, and the resulting value of u\ l had an uncertainty of greater than 100%. 
Similarly, allowing one centrifugal BOB parameter t\ l to vary led to reductions in dd of less that 
0.05%. As a result, the model actually used to obtain the results summarized in Fig. 5 included 
only the one free BOB parameter Uq 1 . 

All else being equal, the "best" model for a given system is the one which achieves an optimum 
quality of fit (lowest dd) with the smallest number of free parameters. When more than that 
minimum number of parameters are used, the additional degrees of freedom in parameter space 
will not be strongly constrained by the data, and the possibility of problems in the extrapolation 
regions tends to increase. On this basis we choose the M3LR|'3(3) model corresponding to the 
'gradiant-shaped' black points at r ref = 8.0 A on Fig. 5 as our preferred model. To four decimal 
places, increasing N by one had no effect on the associated value of dd = 0.7069. While the 
analysis described above led to our determination of optimal models for the a 3 £+-state potential 
and BOB function, the associated parameters were also allowed to vary freely in the global two-state 
fit which simultaneously determined potential energy functions for the 1 3 E+ state. Hence, they 
will be reported later. 

One final component of this discussion is an illustration of the remarkably robust extrapolation 
properties of the MLR potential form. In the initial stages of this study, only the optical 1 3 S+ — 
a 3 E+ data which spanned vibrational levels v(a) = — 7 were considered in the analysis. The 
highest of the associated a 3 £+-state levels is bound by 26 cm -1 for 7 ' 7 Li2 and by 16 cm -1 for 6 ' 6 Li2. 
Nonetheless, the optimal MLR potential obtained from that analysis was an M3LR|*3(r) function 
whose well depth of D e = 333.79(1) cm -1 is very close to the value 333.76(1) cm -1 yielded by 
the analysis of the full a 3 £+-state data set (see Fig. 5). Moreover, the v(a) = 10 binding energy 
predicted by that potential was 0.4222 cm -1 , which is remarkably close to the measured PAS valued 
of 0.4160(±0.0013) cm -1 . Thus, a DPF analysis using an MLR potential with a good multi-term 
theoretical mlrO") seems capable of yielding quite reliable extrapolations to predict both the distance 
from the highest observed level to dissociation and the number and energies of unobserved higher 
levels. 
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C. Model for Li 2 (l 3 £+) and Results of the Two-State Analysis 



Following the discussion of § II.C, the potential energy function for the 1 3 S+ state of Li 2 was 
represented by an MLR potential whose long-range tail was defined by u^(r) of Eq. f )26|) . As usual, 
fits were performed using models corresponding to a variety of values of the exponent polynomial 
order N, of the power q defining the polynomial expansion variable, and of the reference distance 
r ref of Eq. (ITS]) . Since the inverse-power terms contributing to u^(r) have powers m, = {3,6,8}, 
the power p defining the radial variables y^ir) and y T p ei (r) of Eqs. jSJ) and (TT5]) was fixed at p = 
6 (> 8 - 3), and while Cf and Cf were held fixed at the theoretical values of Tang et al.,— Cf 
was treated as a free parameter. All of these fits treated the full range of data for the 1 3 S+- 
a 3 E+ system, and while the potential for the a 3 E+ state was represented by the M3LR|g(r) model 
described in § III.B, its parameters were also free variables in these fits. 

Figure 6 summarizes results for six families of 1 3 £+-state PECs: those for exponent polynomial 
orders N = 6 — 9 with q = 3 being represented by solid points, while those for polynomials orders 
N = 8 — 9 with q = 4 (bottom panel only) are shown as open points. As is expected, the fact 
that the 395 observed 1 3 £+-state level energies are now constrained to be eigenvalues of a potential 
function rather than being free fitting parameters means that the dd values associated with the best 
of these fits are somewhat larger than those for the a 3 £+-state analyses summarized in Fig. 5. 

As has been the case in other treatments of this type, for any reasonable values of q and r re f 
the fits converge to essentially the same optimum dd value when the polynomial order N becomes 
sufficiently large.— i2i>22ilZd2, A manual optimization of r re f is undertaken in order to determine a 
'best' model, which is defined as one which: (i) gives a good fit to all data, (ii) is defined by 
the smallest number of free parameters, and (iii) has no unphysical behavior in the extrapolation 
regions. As is usually the case, models with larger q values (here q , =4) require a higher-order 
polynomial to achieve a given quality of fit. While not shown, the dd values for N — 10, q = 4 
models with r re f = 3.4 — 3.6 are essentially identical to those for N = 9, q = 3 models with the 
same r re f (solid round points), but are bigger at larger and smaller r re f values. Models with q = 1 
or 2 tended to have inflection points on the short-range repulsive wall, even for cases with fairly 
large r re f values. 

The results in the two upper panels of Fig. 6 show that the values of D e and Cf yielded by fits 
to models with N = 9, q = 3 vary relatively slowly with r re f, and that for r re f values which give 
small dd values, other types of models yield very similar results. For both of these properties the 
analogous results for q = 4 models were much more strongly model-dependent. We therefore chose 
the N = 9, q = 3 potential with r re f = 3.6 A as our recommended model for this state. 

The error bars shown in Fig. 6 are the 95% confidence limit uncertainties in the parameter yielded 
by the non-linear least-squares fits. Although the binding energies of the highest observed levels for 
the 1 3 S+ state are even smaller than was the case for the A 1 !]^ state,— the uncertainties in the 
fitted Cf values obtained here (> 75 cm -1 A 3 ) are an order of magnitude larger than the analogous 
uncertainties yielded by the A — X analysis of Ref.llOj. It is not clear why this should be the case, 
other than the fact that the data gap from t>(l 3 £+) = 8 to 61 for 7,7 Li 2 or to 55 for 6 ' 6 Li 2 may 
be expected to introduce additional uncertainty into the analysis of the limiting near-dissociation 
behavior. However, the difference between the Cf value implied by the present analysis and that 
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determined from the A — X analysis of Ref.llO (dash-dot-dot line in the uppermost panel of Fig. 6) is 
significantly larger than the mutual uncertainties. Moreover, repeating the present analysis with Cf 
fixed at the value yielded by the A — X analysis (357829(±8) cm -1 A 3 ) increased the overall value 
of dd by 0.8%, and increased dd for the PAS data by a massive 21%! It may be that a combined 
5-state analysis of the data sets for the two cases will resolve this discrepancy, but that is beyond 
the scope of the present work. Thus, our recommended model for the 1 3 £+ state is an M3LR|J(9) 
potential with mlr(^) defined by Eq. (126j) . and with Cf determined from the fit. 

The parameters defining our recommended models for the a 3 E+ and 1 3 £+ states of Li2 are listed 
in Table II. The fact that the uncertainties in the values of D e and r e are an order of magnitude 
larger for the 1 3 £+ state is to be expected, both because of the 5200 cm -1 data gap for v > 7, 
and because the fact that its lowest observed level is v — 1 means that the extrapolation to the 
potential minimum is much longer for this case. As for the a 3 E+ state, obtaining a good combined- 
isotopologue fit required the introduction of BOB corrections in the effective adiabatic potential 
function for the 1 3 E+ state. As shown in Table II, our recommended model includes a third-order 
polynomial expression for the 'adiabatic' correction radial strength function S^\(r). Increasing this 
polynomial order further or allowing for a non-zero centrifugal BOB function yielded no significant 
improvement in the quality of fit, while reducing this polynomial order by one or two terms increased 
the dd value for the fit by 1.6% or 3.8%, respectively. 



D. BOB Functions and Isotope Effects 

The radial strength functions defining the effective adiabatic BOB correction to the potential 
energy functions for both states are shown in Fig. 7. Since 7,7 Li2 was chosen as the reference 
isotopologue, Eq. flMJ shows that the fact that g^ 3 ^}^ 3 ^} ')- S^ i] (oo) = u { a ^ ±} -u { ^ ±} = 
Mq" is positive means that T)\ a is (slightly) larger for 6,6 Li2 than for 7,7 Li2. The upper curve 
in Fig. 7 shows that the (r) contribution to the isotopologue dependence of De is also 

positive (and much larger). However, the isotopologue dependence of the 2B^ a \r) contribution to 
the effective adiabatic potential (see Eq. ( 1361) ) makes a negative (—0.171 cm -1 ) contribution to the 
difference Se° t(6,6) — S)* ^ 7 ' 7 ) for the 1 3 E+ state, and turns out to be the dominant term. 

The results presented in Table II were obtained from an analysis which treated 7 ' 7 Li 2 as the ref- 
erence isotopologue, and the first row of Table III presents characteristic properties of the resulting 
a 3 E+ and 1 3 £+ potential energy functions for that species. Note that the values of 3)g 0t (l 3 S+) and 
7"g 0t (l 3 S+) were obtained after combining the MLR potential with the additive adiabatic correction 
term of Eq. (1351) . The next two rows of this table then show, respectively, the isotopic changes in 
and the resulting values of these quantities for 6,6 Li 2 , as implied by the BOB correction functions 
Sad(?") and (for the 1 3 £+ state) the 2 B^ a \r) term (see Eqs. (l3TI)-(l36l)). Of course it is equally 
feasible to perform the overall analysis using 6,6 Li 2 as the reference isotopologue, and the last row 
of Table III shows the properties of that isotopologue obtained in that more direct manner. It is 
reassuring to see that within the uncertainties, the results in the last two rows of this table agree 
with one another. 

Of course it is simpler to work with potential functions that do not require the addition of separate 
adiabatic correction functions AV^\r). Hence, for the convenience of those interested primarily 
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in the minor isotopologue 6,6 Li2, a version of Table II for the case in which this species was used 
as the reference isotopologue is included in the Supplementary Data supplied to the journal's www 
archive. 



IV. DISCUSSION AND CONCLUSIONS 

A combined-isotopologue DPF analysis of 2692 optical data for the 1 3 £+ — a 3 X+ and 2 3 II 9 — 
a 3 E+ band systems of 7,7 Li 2 and 6 ' 6 Li 2 , together with 99 PAS data for the 1 3 S+ state and one for the 
a 3 S+ state, has yielded analytic potential energy functions for the 1 3 E+ and a 3 S+ electronic states 
which (on average) explain all of those data within the experimental uncertainties (dd = 0.789). 
The present potential energy function for the a 3 E+ state of Li 2 is one of the most accurate ground- 
triplet-state potentials determined for any alkali-atom pair. The resulting scattering lengths for 
7,7 Li 2 and 6 ' 6 Li 2 are osl = —14.759(9) A and —1906(50) A, respectively, where the uncertainties 
were estimated by repeating the overall analysis with C 6 3 increased/decreased by 0.01% from 
the recommended values of Tang et alM This 0.01% is a factor of 3 larger than the Cg uncertainty 



reported in Ref.|42|. The uncertainty in osl is much larger for 6,6 Li 2 than for 7,7 Li 2 simply because 
all else being equal, scattering lengths that are very large in magnitude are much more sensitive to 
details of the potential energy function. Listings of band constants (G v , B v , D v , H v , etc.) calculated 
from this potential for all bound levels of the a 3 E+ state for all three Li 2 isotopologues are included 
with the Supplementary Data supplied to the journal's www archive. 

For the 1 3 E+ state, the present analysis has provided an analytic potential energy function which 
very robustly bridges the ~ 5200 cm -1 gap shown in Fig. 1, between the fluorescence measurement 
domain v = 1 — 7 and the PAS data for v > 56 for 6,6 Li 2 and > 62 for 7,7 Li 2 . To illustrate this 
interpolation behavior, Fig. 8 plots calculated properties of our recommended potential for 7 ' 7 Li 2 
in the manner suggested by near-dissociation theory (NDT).— In particular, NDT predicts 
that for vibrational levels lying near the dissociation limit of a potential whose limiting long-range 
behavior is defined by an attractive C3/7" 3 interaction energy, the l /& power of the binding energy 
(D — E v ), the l /§ power of the vibrational level spacing AG„ + i/ 2 , and the x /a power of the inertial 
rotational constant B v , should all be linear functions of v , with slopes determined by the value of 
the C3 coefficient. The solid triangular points in Fig. 8 represent the experimental data, while the 
open round points are our predictions for the 'no-data' regions. The dash-dot-dot lines in Fig. 6 
are the limiting NDT slopes implied by the fitted Cf value of Table II. The deviation from this 
behavior at very high v reflects the fact that the 3x3 interstate coupling reduces the magnitude 
of the effective C3 coefficient in the limiting region by a factor of 1 /3 (see Fig. 4) as one approaches 
the limit. Calculated band constants for all bound levels of all three Li 2 isotopologues in this state 
have been placed in the journal's Supplementary Data archive. 

It is noteworthy that predictions generated from a variety of other MLR potential models (i.e., 
models defined by different N or q values) which yield good fits to the data are identical on the 
scale of Fig. 8. This model-independent bridging of a data-gap spanning 73% of the well depth 
is a remarkable illustration of the robustness of the MLR potential function form. The ability of 
this function to readily incorporate the effect of two-state^ or three-state (present work) coupling 
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in the long-range region is a further demonstration of its capabilities. A Fortran subroutine for 
generating the recommended potentials is one of the items placed in the journal's Supplementary 
Data archive. 

One puzzle left by this work is the discrepancy between the value of Cf for interacting Li( 2 P) + 
Li( 2 S*) atoms determined in the present analysis (3.575 57(78) x 10 5 cm -1 A 3 ) and those obtained 
in the v4( 1 S+)-state analysis of Ref.Eol (3.578 29(7) x 10 5 cm -1 A 3 ) or from the recent theoretical 
calculations of Tang et al.& (3.578108 9(7) x 10 5 cm^A 3 ). While small on an absolute scale, 
this 0.076% discrepancy is much larger than the estimated uncertainties, and repeating our overall 



analysis with Cf fixed at the A-state value from Ref.llOl yielded a distinctly poorer quality fit, 



especially for the PAS data. It may be that a combined five-state analysis of all of the data 
considered here with those used in the A — X analysis of Ref.llOl will shed light on this question, 
but that will have to await future work. 
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Table I: Summary of experimental data used in the present work. 



isotop. 


type 


unc. (cm 1 ) 


v(2 3 U g ) 




%+) 


v{a^+) 


* data 


source 


7 < 7 Li 2 


1 3 £+ emission 


0.01 


— 


1 


-7 


0-7 


1279 


Ref. .12 




2 3 II 9 emission 


0.005 - 0.01 


1-2 






0-9 


137 


Ref. 14 




PAS(1 3 £+) 


0.0043 — 0.00073 




63 


— 90 




30 


Ret. 16 




PAS(a 3 £+) 


0.0013 








10 


1 


Ref. 15 


6 ' 6 Li 2 


1 3 £^~ emission 


0.01 




1 


- 7 


0-7 


1276 


Refs.13 




PAS(1 3 S+) 


0.00110 - 0.01067 




56 


-84 




69 


Ref. 16 




Overall 




1 - 2 


1 ■ 


-90 


0-10 


2792 
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Table II: Parameters defining the recommended MLR potentials and BOB functions for the a 3 £+ 
and 1 3 £+ states of Li 2 obtained using 7,7 Li 2 as the reference isotopologue. Parameters in 
square brackets were held fixed in the fit, while numbers in round brackets are 95% confidence 
limit uncertainties in units of the last digits show. The analysis used the 7 Li 2 Pi/2 <— 2S , 1 / 2 
excitation energy of 14903.648130 cm -1 and 2 P%/2 ^— 2Pi/2 spin-orbit splitting energy of 
0.335338 cm -1 from Ref.l38l Units of length and energy are A and cm -1 ; the exponent 
expansion coefficients 0i are dimensionless, while the parameters U( defining the 'adiabatic' 
BOB strength function of Eq. (12 7p have units cm -1 . 





«( 3 ^) 




c(l 3 £+) 




333 758 (7) 




7f)Q3 44 f3l 


T 

1 e 


4 17005 (3) 




3 06514 (91 


Cc 

^ 


[6 7185 x 10 6 1 




3 57557 (78") x 10 5 


^8 


\1 12629 x 10 s 1 


^6 


[1 00054 x in 7 1 




\2 78683 x 10 9 1 


^8 


[3 69953 x 10 s 1 


Dt ■ 

rLi 


r 541 




N/A 


{P, <?} 


{5, 3} 




{6, 3} 


r re f 


[8.0] 




[3.6] 


Po 


-0.51608 




-1.6373863 


01 


-0.09782 




0.29197 


02 


0.1137 




-0.55544 


03 


-0.0248 




-0.2794 


04 


— 




-1.5993 








—0 673 


06 






-1.23 


07 






-1.29 


08 






0.5 


09 






2.6 


{Pad, <?ad} 


{6, 6} 




{3, 3} 


U 


0.059(11) 




1.367(7) 


Ml 






2.7(4) 


U 2 






-1.3(5) 


U 3 






-1.8(4) 


Uoo 


[0.0] 




[1.055740] 
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Table II. B: (for Supplementary Data submission) 

Parameters defining the recommended MLR potentials and BOB functions for the a 3 £~ 



and 



1 3 S+ states of IA2 obtained using 6,6 Li2 as the reference isotopologue. Parameters in square 
brackets were held fixed in the fits, while numbers in round brackets are the 95% confidence 
limit uncertainties in units of the last digits show. The analysis used the 6 Li 2 P\/2 ^— 2S1/2 
excitation energy of 14903.296792 cm^ 1 and 2 P 3 /2 ^~ 2P1/2 spin-orbit splitting energy of 



0.335324 cm -1 from Ref.l38l Units of length and energy are A and cm -1 ; the exponent 
expansion coefficients B% are dimensionless, while the parameters Ui defining the 'adiabatic' 
BOB strength function of Eq. (12 7|) have units cm -1 . 





a( 3 K) 

V U / 




c(l 3 E+) 




333.778 (7) 




7093.54 (3) 


r e 


4.17001 (3) 




3.06520 (9) 




[6.7190 x 10 6 ] 


c 3 


3.57549 (78) x 10 5 


C s 


[1.12634 x 10 s ] 


c? 


[l. 00059 x 10 7 ] 


wo 


[2.78694 x 10 y J 




[3.69965 x 10 s J 


Phi 






NT/A 

IN / i\ 


{"Pi q\ 


{5, 3} 




{6, 3} 


r re f 


[8.01 

L J 




[3.61 

L J 


Po 

1 u 


-0.516086 




-1.6373493 


61 

1 -L 


-0.09783 




0.29191 


02 


0.1136 




-0.55591 




-0.0249 




-0.2788 








-1.5931 


fa 






-0.685 


Be 






-1.306 


07 






-1.34 


08 






0.6 


09 






2.7 


{Pad, gad} 


{6, 6} 




{3, 3} 


U 


0.067(13) 




1.595 (8) 


Ml 






3.1(5) 


"2 






-1.5(5) 


u 3 






-2.1(5) 




[0.0] 




[1.23141] 
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Table III: Properties of the recommended potential energy functions for the a 3 S+ and 1 3 S+ states of 
Li 2 (energies in cm -1 and lengths in A), with 'changes' calculated using Eqs. (|3ip - (|36p . The first 
three rows correspond to use of ' Li 2 as the reference isotopologue, while for the last row it was 
6 < 6 Li 2 . 





a 3 S+ 

U 


state 


AT e 


1 3 £+ state 


fit isot. 




r e 


/j"\tot ^tot 


2-isot 7 ' 7 Li 2 


333.758(7) 


4.17005(3) 


8144.989(43) 


7092.417(33) 3.06524(9) 


change 


0.019(4) 


-0.000015(3) 


-0.265(5) 


-0.067(2) 0.000075(1) 


6 < 6 Li 2 


333.777(8) 


4.17003 5 (3) 


8144.724(43) 


7092.349(33) 3.06532(9) 


2-isot 6 ' 6 Li 2 


333.778(7) 


4.17001(3) 


8144.726(33) 


7092.347(33) 3.06532(9) 
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Figure 1.: Overview of the potential energy functions and observed vibrational levels of the l 3 Sjj~ and 
a 3 S+ states associated with the present analysis. The insert shows a fragment of the 1 3 S+ state 
potential at the energy range associated with the 6 ' 6 Li2 PAS data. 
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Figure 2.: Eigenvalues of Eq. (|18p for the l g states of L12 dissociating to the Li(2s 2 5') +Li(2p 2 P) limits, 
as generated using the long-range coefficients of Tang et air- 



25 




26 




Figure 4: Comparison of four representations of the long-range potential for the l 3 Sj~ state of Li2, with 
energies in cm -1 and lengths in A. 
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Figure 5: Determination of the optimum model for the a 3 X+ state of L12. 
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Figure 6: Determination of the optimum model for the l 3 Sjj~ state of Li2. 
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Figure 8: Spectroscopic properties of the 1 3 S+ state of 7 ' 7 Li2. 



